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Abstract 

An isospin and spin dependent form of the equation 
of state for nuclear matter is presented. This form 
is used for the description of nucleon interaction in a 
new dynamic model. Preliminary calculations show 
that the new approach makes possible predicting the 
alpha structures appearing in the case of the ground 
state even-even nuclei. The change in the structure of 
nuclei built on the basis of alpha particles occurring 
for ^°Zn is discussed. 

1 Introduction 

Although nuclei are close to 55 orders of magnitude 
lighter than a typical neutron star, it seems that the 
state of the material comprising of the above two ob- 
jects can be determined by the same equation, which 
is the equation of state (EOS) of nuclear matter. In 
this case EOS is defined as the average energy per 
baryon expressed as a function of thermodynamic 
variables. It would seem that the above statement 
can be more easily justified for heavy nuclei but, as 
can be shown, using an appropriate form of the equa- 
tion of state can give good results also for the descrip- 
tion of the ground states for lighter nuclei. We arrive 
at these conclusions using semiclasical, microscopic 
version of the liquid drop model in which dynamics 
is governed by EOS. 

The description of dynamics of nuclear system, 
which uses the concept of the equation of state, can 
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be found in the reaction models referring to the liq- 
uid drop model (LDM) \1^. In this description, it is 
assumed that the energy of the system is determined 
mainly by using the following three components: 

- the volume energy defined by the kinetic energy 
related to internal, fermionic motion and respective 
potential interactions for infinite nuclear matter in 
equilibrium, 

- the surface energy, which in classical description 
could be treated as a result of the surface tension 
action, 

- the Coulomb energy associated with the proton 
charge. 

Such models typically assume one body dissipation 
and no compression of the nuclear matter. In such 
an aproach the volume energy together with energies 
which are corrections related to a not equal number 
of protons and neutrons and a non-zero spin of the 
system may be treated as the 0*'* term of the expan- 
sion of the energy density function (given by EOS). 
No matter compression assumed in this model intro- 
duces additional constraints on the global coordinates 
which describe the system shape and its evolution. 

In our microscopic description of the system evo- 
lution, we take into account also the next term of 
the expansion of the equation of state as a function 
of the density of nuclear matter. The presented ap- 
proach treats the local density of nuclear matter as a 
result of the aggregation of partial densities given by 
Gaussian packets which represent the nucleons form- 
ing the system. As a result, the local density of mat- 
ter is explicitly given by the position of these packages 



and their variance. In our approach we describe the 
volume energy as the respective functional (volume 
integral of the energy density given by the equation 
of state multiplied by the density of the nuclear mat- 
ter). 

The nuclear EOS describes equilibrated, infinite 
nuclear matter and is necessary to calaculate the 
properties of astronomical objects as, for example, 
neutron stars or supernovae [2]. In such large objects, 
the surface effect in the description of the total energy 
is negligible. For small objects, such as atomic nuclei, 
the energy due to interactions at the surface, cannot 
be ignored and in such case the Hamiltonian values 
must have additional contributions, corresponding to 
this additional surface energy. In the presented ap- 
proach, the surface energy is related to the dispersion 
of the energy density in the nucleon's environment 
and also to the increase of the diffuseness of the nu- 
cleon momentum caused by nucleon location on the 
surface (reduce the uncertainty given by surface in- 
teractions as described in [3]). 

In the interaction description which is applied in 
the model presented here one can distinguish ingre- 
dients related to the volume, surface and Coulomb 
energy, which are the essence of description used in 
LDM model . 

In our approach, all the aforementioned types of 
energy are functionals expressed by integrals deter- 
mined by the distribution of nuclear matter. Since 
these densities are here uniquely determined by the 
position and other parameters of wave packets (spin, 
isospin and variance), the description can be inter- 
preted as the Microscopic LDM model (MLDM). 

As we know, the equation of state of nuclear matter 
is still not well established and its determination is 
one of the most important tasks faced by experimen- 
tal and theoretical nuclear physics. In the presented 
approach, the equation of state plays a very impor- 
tant role and in this way the described model has an 
opportunity to contribute in EOS investigation. It 
is also important to use a convenient form of EOS 
suitable for efficient description of the heavy ion re- 
action as well as for description of the ground state 
properties of the reacted nucleii. This form of EOS 
should fulfill the following requirements: 

1. The compression part of the equation of state. 



in the neighborhood of saturation density po (where 
Po is the density of the matter which is equilibrated in 
isospin and spin and for which the density of energy 
reaches a minimum) can be written as the sum of the 
two types of energy: 

- eo: energy related to symmetric nuclear matter 
SNM; 

- esym- energy produced by the nuclear matter 
symmetry disorders, which is called symmetry energy 
(see section 2). 

It turns out that by the adoption of certain sim- 
plifying assumptions (see appendix A), these energies 
can be expressed as a function of the total density of p 
. The expectation which is in line with a mean field 
theory (e.g. with the Skyrme-Hartree-Fock model) 
determine the zero energy density in case where den- 
sity of matter tends to zero (it is in contrast with 
some other models like, for instance, an idealized 
alpha-mater picture [1]). 

2. For the saturation density we require: 

- specific value of the energy density; 

- specific values of the first and the second deriva- 
tive of the energy density (defined by EOS); 

and for density p ^)- the additionall requirement 
is e{p) -> 

then the simplest form of a polynomial, which can 
meet all these requirements can be expressed as: 



^ap + 13 p^ + 7p^ 



(1) 



The method of determining the coefficients alpha, 
beta, gamma is given in Section 2. and Appendix 
A. As it is shown [3], form ^ is suitable for differ- 
ent theoretical descriptions of EOS. It is also worth 
noting that this form of the EOS allows us to find 
the respective mean energy values in an analytical 
way (in case when the density is expressed by partial 
density, taken as the Gaussian functions). 

A modified equation of state in a form suitable for 
MLDM calculations is defined in the next section. In 
the third section we will present the application of 
the proposed form of the EOS together with an ad- 
ditional surface energy term used in MLDM model. 
The obtained preliminary results using this model are 
related to the properties of nuclei in its ground states 
and are presented in Section 4. Since in our descrip- 



tion we use EOS, in which energy symmetry is as- 
sociated with spin and isospin asymmetry (see next 
section), the geometric structure defined by the po- 
sition of the centers of the wave packets can be in- 
terpreted as clusters. For the nuclei with even and 
equal number of neutrons and protons these clusters 
appear as alpha particles. In Section 5 we present 
such alpha structures, which appear in model calcu- 
lations for light nuclei. 

In the last chapter conclusions and future plans are 
presented. 



2 An equation of state given by 
third-degree polynomials 

The concept of the equation of state of nuclear matter 
in the case of two-component proton-neutron gas has 
been discussed in many articles. This equation as well 
as a form based on the use of the description of third- 
degree polynomials is described in [3] and works cited 
therein. In the next section we present the ground 
state properties of nuclei which are obtained by the 
microscopic, density-functional approach. EOS plays 
a key role in this approach. 

This model is designed to explore the dependence 
of the dynamics on wave packets, isospin and the 
spin of nucleons, which these packets represent. By 
isospin dependence, we mean the dependence on the 
charge of the nucleon (third component of isospin vec- 
tor). 

If the isospins, spins and variances of wave packets 
are the variables affecting the mutual interaction of 
nucleons, then appropriate equations of state must 
describe a gas consisting of four components: 

-protons with spin up - density pp^ 

-protons with spin down - density ppj, 

-neutron with spin up - density p„-|- 

-neutron with spin down, density pn^ 

In our considerations, by 'nucleons with spin up' or 
'nucleons with spin down' we mean nucleons having 
spin projections 'up' or 'down' on the chosen quanti- 
zation axis. Now the EOS can be formally written: 



According to symmetries characterizing the nuclear 
interaction, the system energy is conserved during 
mutual exchange of neutrons and protons, and also 
when changing the projection of the spin of all par- 
ticles to the opposite spin. To allow us to use these 
symmetries we define new coordinates: 
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where p is the total nuclear matter density, and po 
is the density of isospin and spin balanced matter at 
saturation. We find that during the above-mentioned 
operation (mutual conversion of neutrons into pro- 
tons and the change all the nucleon spin projec- 
tions) the sign of the coordinates: neutron-proton 
assymetry 5 and spin assymetry ?7„ and rjp are re- 
spectively changed. Therefore in order to assure en- 
ergy conservation, in the expression which defines the 
system energy the values 5 and rjn and rjp can occur 
only as products with the respective even number of 
factors. Using these symmetries one can show (see 
Appendix A) that the equation of state can be writ- 
ten with the components having an extended form of 
symmetry energy; 

e = eoo + —4 + 



5^ 



e/0 + -TT^ 



18^ 



ril) 






2?7„ryp Cijo 



'<- 



K., 



18 ^ 



(7) 



e(Ppt'Pp4.'Pnt'Pni) 



(2) 



The first two components of the sum ([7|) represent 
a well-known form of the equation of state. The first 
component describes the matter in a balanced system 



(zero isospin and spin), where Kq is the coefficient 
of compressibility of nuclear matter and the second 
one describes the isospin interaction is given in the 
following form: 



6/ = e/o + -—i 



18 



e 



(8) 



which is called the symmetry energy for which e/o is 
the Wigner constant and the coefficients L/ and Ki 
are the slope and curvature respectively. 

In equation ([7]) , by analogy to the isospin symme- 
try energy, one can distinguish the spin symmetry 
energies for the neutrons and protons separately: 



and energy: 






e. = e,o + ^e + ^e 



(9) 



(10) 



for the mutual, spin interaction of protons and neu- 
trons. Factors develop spin symmetry energy: euo, 
La and Ku, (constant slope and curvature) with in- 
dices ii, describe the energy relating to the proton 
or neutron gas. While the respective indices ij indi- 
cate that the symmetry energy refers to the mutual 
interaction of protons and neutrons. 

All types of the symmetry energy describe the im- 
pact of the spin and isospin polarization of the matter 
on the average energy of the nucleon (EOS). Accord- 
ing to the observations presented in the introduction 
of work, we assume that the energy associated with 
the SNM as well as all types of symmetry energy 
should disappear as the density of matter tends to 
zero. Therefore, to ensure zero values for all energies 
occurring in equation ([7]) at p = and to get the 
proper slope and curvature of the saturation density 
Po we use the form of third-degree polynomial pro- 
posed in [5] . As mentioned in the cited work, if the 
density of matter is described by the sum of Gaussian 
distributions, then such a form of the EOS allows us 
to calculate the average energy analytically. 

Any type of energy found in ([7]) may be written as: 



and the index 'k' can take values 0, /, ii, and ij. 
Here the coefficients a^, /3fc, jk are determined by 
the relations: 

ak = f 3efco - -Lk + T^j /Po (12) 

3efco -Lk + ^) IpI (13) 

IpI (14) 



lk = efeo 



3 18 



3 Application of the proposed 
EOS form to the MLDM cal- 
culations 

The proposed MLDM model, similarly like in 
M 1 Hn] : [H] : describes the time evolution of the wave 
function represented by the product of M Gaussian 
wave packets, which represent the nucleons forming 
the system : 



M 



$=n' 



(15) 



fe=i 



i2na] 



2^3/4 



exp 



ir-{rk)r 
^4 



-r{Pk) 



(16) 

where af, (vk), {Pk), are the width of a Gaussian 
wave function and the centers of its position and as- 
sociated momenum for each of M nucleons. In our 
approach the variables (rfc), (p^) and a'^ are time- 
dependent parameters describing the wave functions. 
Additionaly, we assume that with every wave packet 
for each of M nucleons, isospin and spin are associ- 
ated and that they remain fixed during the interac- 
tion. 

The equations of motion of variables (r^), (pk) and 
al are derived using the time-dependent variational 
principle (see e.g. [S]) based on the action minimiza- 
tion. For this purpose, we define the action as: 



S = 



e/c = akP + l3kP^ + IkP^ 



(11) 



t2 

I L{^,^*)dt (17) 



with the Lagrange functional given as: 



dt 



dt 



L = ( $ ifi— - H $ ) = ( $ ifi— $ ) - ($ \H\ $) 

(18) 

Now we will concentrate on the determination of 
average value of the Hamiltonian($ \H\^). As usual 
the Hamiltonian can be written as the sum of the ki- 
netic energy T and potential interaction V . To deter- 
mine the average kinetic energy we assume that each 
wave packet has a momentum dispersion Up^^ around 
its mean value (pfc). Then the average kinetic energy 
associated with the system: 



k=M 

fe=l 
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2m 



(19) 



We assume also that the interaction between nucleons 
can be described by the potential given in form; 

V ({rk}) = Vn ({rk}) + Vs ({rk}) + Vco^a (K}) 

(20) 
where Vn describes the nuclear interaction, Vs de- 
scribes the modification of the interaction induced 
by changing the density of matter around the nucleon 
(surface energy) and Vcoui consists of the Coulomb 
interaction. If the nuclear matter is characterized 
by non zero spin polarization then, in the space a 
magnetic field can occur. This field can also be gen- 
erated by the movement of protons. In such a case 
the Hamiltonian should be corrected for the magnetic 
interaction. 

Now the averege value of the Hamiltonian can be 
expressed by: 



{^\Vn\^) + {^\Vs\^) 



k=M 
(<I>|H|$)= ^ 

{<^\Vc. 



(PkT 

2m 



^< 



2m 

($|Af|<I>) 
(21) 



In our further considerations we ignore magnetic in- 
teraction. The crucial assumption of the present ap- 
proach is that energy of the fermionic internal motion 



can be described by 



3<7i 



2m 



SO the sum 



Bv = 



k=M 

E 

k=l 



2m 



{^\Vn\^)= / e(p,^,cr„,crp)p(r)d'^r 

(22) 

can be treated as a volume term of binding energy. As 
we know, the ground state energy of nuclear matter 
is described by the EOS, which considers the energy 
density as the sum of kinetic energy, associated with 
the internal fermionic motion, and energy given by 
the potential interactions. Therefore in our approach 
we do not have to deal with finding the appropriate 
distribution of the nucleon momenta, provided that a 
correct description is given by the selection of appro- 
priate EOS parameters. Such an EOS parameteriza- 
tion replaces the potential parameterization normally 
used in such approaches. 

As in ||3:, the surface energy is defined as 



B 



surf 



{^\Vs\<^)^so 
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(23) 



where in this semiempirical part sq is the surface en- 
ergy coefficient. In this formula, the variance (Te(fc) 
denotes the variance of energy e{p, S, rjnjVp) with re- 
spect to the probability ^p{r) of finding the fc-th nu- 
cleon: 



(Je{k) 



2 k„ 



{e~e{p,S,T]n,Vp)) ■pir)d'^r 



(24) 



where e is the everage energy given by the formula: 
e(fc) = / e(p, S, 7?„, 77p) ■''p{r)d^r (25) 



In appendix 1 we describe the metod of calculation 
used for variance ae{k) calculations. 

Let us discuss briefly the description of the sur- 
face energy given by formula (P^. In our approach 
we consider drops of nuclear matter shown in Figure 
1. The density of matter p{r) — /9„(r) + Pp{r) is a 
function of the distance r and pn (r) — pnf + Pn\. and 
Pp('') = Ppt + Ppi- If ^^'^ ^*'' packet is in the region 
of fixed energy e(p, 5, r]n, ijp), then its variance cre(l), 
associated with this package ''p{r), is equal to zero. If 
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Figure 1: An example of the position of the partial 
density function and associated energy variance. As 
one can see for a position close to the center, variance 
of the energy ((24|l vanishes. See text for more details. 

the package V(i") is in the change of the energy, then 
the variance Ueik) is greater than zero. According to 
formula (|23|) . the surface energy associated with the 
i*'' packet is inversely proportional to the uncertainty 
of location, which means that it is proportional to the 
indeterminacy of the package momentum. Now for- 
mula p3|) can be understood as follows. Variation 
of energy ae{i) (and associated forces) describes the 
possibility of localization of the i*''-th nucleon in case 
of its entering the alternating field. Such localiza- 
tion forces result in an increase of indeterminacy of 
the momentum and an increase of the average kinetic 
energy associated with the i*'* package . 

The last term in (|20p is the Coulomb energy (for 
protons only) and is given in the form: 

($ \Vcoui\^) - \Y. hp^^)r^\ W)d'vd'r' 
2^y |r-r'| 

(26) 
This integral can be calculated analytically as (see 
eg. [6]): 



where rki = |(rfc) - (r;)| and sigma a = ^J a\ + a'f . 

Finally, the expected form of the Hamiltonian can 
be written in a form which is analogous to the de- 
scription given by the liquid drop model (LDM) : 

(* 1^1 *) = E %^ + By ({(^^)} ' {-^» + .28) 

Bs ({(rfe)} , K}) + Vcoui {{{vk)} , Wk]) 

Based on the Euler-Lagrange equations one can 
find the corresponding equations of motion for co- 
ordinates 

(Pfe) , (r/c) and a^ [see Euler-Lagrange equations]. 
This will be described in a forthcoming paper. 

Here, as the first test of the new form of the equa- 
tion of state, based on the MLDM model, we will try 
to describe the ground state properties of a selected 
set of nuclei. In our very preliminary approach we 
will focus on binding energiess and nuclear charge 
radii (RMS). 

4 The MLDM calculations for 
describing the ground state 
properties of nuclei 

In the MLDM model the ground state of a nucleus 
is a many-body state which is an absolute minimum 
with respect to variation of (r^) , ak- Additionally, in 
the ground state, for every k the average momentum 
{Pk) — , which means that for MLDM the particle 
will be essentially at rest and the system corresponds 
to a solid. In order to determine the nuclear ground 
state wave function the following procedure of search- 
ing the set {(rfc) ,afc} is used: 

1. In a limited space we choose a random posi- 
tion of the wave packet together with a certain initial 
variance. In order to accelerate the calculations, for 
even-even nuclei, we can also assume a correlation 
(mating) for packets describing the nucleons of the 



same type and differing only in the direction of the 
spin projection. In such a case we assume that the 
positions of centers and the variances of a chosen pair 
are equal. 

2. Using a selected set of parameters describing 
the interaction (type EOS) we search to find such 
values of these parameters for which the energy of 
the system is minimal. Here, a variety of methods 
may be used. One should also be careful to avoid 
emerging local minima. Wave functions found in this 
way are characterized by certain symmetries, which 
may be helpful in assessing the resulting minimum. 
In the next chapter we will discuss types of symmetry, 
emerging especially for even-even nuclei with an equal 
numbers of protons and neutrons. 

To use the presented model for calculation of 
ground state properties of nuclei, we have to select 
specific values of parameters describing the equation 
of state. For this purpose a global search for optimal 
parameters values should be used. This will be done 
in the next paper, here we apply only a very simpli- 
fied search method based on following assumptions: 

In the first step we try to dctcrmin the p^, ego, Kq 
and So parameters . In our preliminary estimations, 
we assume that po — 0.159 nucl/f?Ti'^ and K — 300 
MeV, which have to be in an experimental and theo- 
retical justification. Then we use the binding energies 
and RMS radii for ■*7Je and ^^C nuclei to determine 
eoo and sq by the trial and error method. In this 
way we obtain ego — —12.9 MeV and sq = 0.09 im? . 
For these nuclei, due to the linking of matter in clus- 
ters which are equilibrated in spin and isospin co- 
ordinates, corrections arising due to the symmetry 
energies appear to be negligibly small. Note that the 
obtained parameters eoo and so, differ from standard 
values (in particular, the typical value of eoo = — 16 
MeV) . This is a result of their assessment for of very 
light nuclei. In the future it will be necessary to de- 
termine these parameters based on a larger database 
(the set of the energy and radii of the nuclei in their 
ground states). 

In the next step, we use the binding energy and 
RMS radii of ^7f , ^7f , ^ife, to obtain an estimate 
of the isospin and spin polarization parameters. Our 
rough estimate gives e/o = 30 MeV, i/ ~ 123 MeV, 
Kj = 500 MeV for isospin parameters and cuq — 



65 MeV, Li, = 300 MeV, Ku = 300 MeV for the 
spin interaction in gas of neutrons (or protons). The 
mutual spin interaction of protons and neutrons are 
described by the parameters e^o = — 1 MeV, Lij = 



20 MeV, K,. 



300 MeV. 



With such parameters, the dependence of the en- 
ergies eo, e/, Cii and e^j, on the density are presented 
in Figure 2. 




n.l 0.2 fl.3 0.1 0.2 0.3 



p [nucleons/fm'] 

Figure 2: The dependence on the density of nuclear 
matter: eo in an isospin and spin balanced state, 
the symmetry energy e/ for isospin of the system, 
the symmetry energy en for the spin polarization for 
proton or neutron matter and the symmetry energy 
energy eij for spin polarization for proton-neutron 
matter (see text) 

In this simple analysis the ^i/, ^i/, ^He, ^_ffe and 
^^C nuclei can be treated as generators of the EOS 
parameters. Reproduced by MLDM model, the bind- 
ing energy and RMS radii are for this generatos ar- 
ranged together with experimental data in Table 1 

For such choice of the EOS parameters (and the 
coefficient of surface tension so — O.OOfm^) we obtain 
predictions for the binding energies per nucleon and 
RMS charge radii these are presented together with 
experimental data in Table 2. 

The data in Table 1 an Table 2 should be trated 
as very preliminary predictions of the MLDM code 



Table 1: MLDM calculation results (for EOS param- 
eters selection, preliminary). Binding energies and 
RMS charge radii of nuclei are compared with the 
experimental data. 



nuclei 


binding 

energy 

data 

(MeV) 


binding 
energy 
MLDM 

(MeV) 


RMS 

charge 

radii data 

(fm) 


RMS 

charge 

radii MLDM 

(fm) 


■'H 


-1.112 


-1.104 


2.14 


2.196 


'H 


-2.827 


-2.809 


1.759 


1.824 


'He 


-2.572 


-2.571 


1.945 


1.974 


^He 


-7.074 


-7.062 


1.676 


1.727 


VJ^ 


-7.68 


-7.719 


2.47 


2.466 



Table 2: MLDM calculation results (prediction, pre- 
liminary) . Binding energies and RMS charge radii of 
nuclei are compared with the experimental data. 



nuclei 


binding 

energy 

data 

(MeV) 


binding 
energy 
MLDM 

(MeV) 


RMS 

charge 

radii data 

(fm) 


RMS 

charge 

radii MLDM 

(fm) 


ibQ 


-7.976 


-7.84 


2.701 


2.613 


^"iVe 


-8.032 


-7.908 


3.005 


2.774 


''Mg 


-8.261 


-7.953 


3.056 


2.874 


3^5 


-8.493 


-7.929 


3.261 


3.127 


^Mr 


-8.520 


-7.896 


3.39 


3.227 


■'^Ca 


-8.551 


-7.85 


3.476 


3.336 



based on a roughly chosen equation of state. For 
heavier nuclei one can see quite large deviations from 
experimental data. In order to achieve greater con- 
sistency with the description of experimental data, 
all EOS parameters need to be re-adjusted. 

5 The alpha structure of even- 
even nuclei with equal num- 
bers of protons and neutrons 

The concept of nuclear clusters appeared along with 
the quantum description of nuclei (and even before 
the discovery of the neutron in 1932 by James Chad- 
wick). In article [7] many experimental arguments 
are presented for the possible emergence of clusters 



of alpha, as a matter substructures in forming nuclei. 
We present here only two of them, both taken from 

i) The binding energy of even-even nuclei with an 
equal number of protons and neutrons appears to be 
a linear function of the number of bonds in alpha- 
particle model, where the number of bonds is equal 
to: fc = 1 for ^Be, fc = 3 for ^^c, fc = 6 to ^^O, for 
fc = 9 ^"Ne, etc. (see Fig. 9). In the cited work this 
phenomenon is explained on the basis of the proposed 
form of the potential interaction between the alpha 
particles, however there is no justification for the in- 
crese in the number of bounds, with an increasing 
number of alpha particles. 

ii) The binding energy of nuclei with an even num- 
ber of protons Z and an odd, higher by 1 number of 
neutrons N = Z +1. In this case we assume that the 
excess neutron is interacting with multi-centrist po- 
tential defined by the system of the appearing alpha 
particles. Let the binding energy for each nuclide X 
be defined by the function 6(X). If the neutron bind- 
ing energy in ^He is denoted by B, then we get the 



following scheme: 
b (^He) 



[b (^Hc) + b (n)] = B 



b (^Be) - [6 f Be) + b (n)] =B + {R + Q) 



b ("C) - [b (i^c) + b (n)] =B + 2{R + Q) 



b (1^0) - [b (i^O) + b (n)] =B + 3{R + Q) (29) 

where the constants R and Q are the result of rea- 
soning, which below is presented in the case of ^Be. 
If a I "01 > and | '02 > denote the neutron wave func- 
tions corresponding to the interaction with the given 
alpha particle, then the two-centre wave function for 
^Be can be expressed approximately as a linear com- 
bination of the single center ones (| ■01 > and | "02 >)■ 
Now, the average binding energy of a neutron, in two- 
cenre potential, can be expressed as: 



b{^Be)-[b{^Be) + 6 (n) 



<^^Jl+i^2\T + Vl + V2 



Fig 3-6 the results of calculations for the nuclei ^^O, 
'^'^Mg, 36 Ar and *°Ca. 

tpl+V2> 



< V"! + ■02 I ■01 + 

(30) 
where in the Hamiltonians H = T + Vi + V2 op- 
erator T describes the kinetic energy and the Vi and 
V2 are related the neutron interaction with the corre- 
sponding alpha particle. As can be seen, the following 
expressions occur twice in the numerator: 

< 01 \ T + Vi I 01 >= B, which is the neutron 
binding energy in ^He, and 

< 01 I V2 I ^01 >= i?, which describes the ad- 
ditional energy associated with the presence of the 
second alpha particle, and 

< V'l \ H \ ip2 >= Q 'AS the energy associated with 
the exchange process. 

If the wave functions | ipi > and | ip2 > sue nor- 
malized, then the denominator in (|30|) is equal to 
2(1-1- < ipi I -02 >)■ Usually, it can be assumed that 
< V'l I 0^2 >= 5" is of little value compared to 1, and 
therefore ([50)1 tends to the corresponding part of the 
expression ((29l) . 

These examples show that the nuclear matter in 
the nuclei has a structure in which the alpha particles 
seem to play an important role. 

The report [7] describes also some of alpha cluster 
models. Our approach is based on single nucleons and 
in this way it is similar to QMD (Quantum Molecu- 
lar Dynamics) [3] , AMD ( Antisimmetrized Molecular 
Dynamics) [TU] and FMD (Fermionic Molecular Dy- 
namics) [llj . in which the dynamics description is 
based on variational principle. 

Since our model includes isospin and spin inter- 
actions occuring between different types of nucleons 
with densities pp-f, ppi, pnf, Pni, nucleons are grouped 
in order to minimize this additional energy. As we 
know, in the case of fermions this causes the connec- 
tion of particles into pairs. For nuclear matter an ad- 
ditional positive energy is asocciated with the lack of 
isospin balance. To minimize both of these energies, 
nucleons create alpha clusters. This phenomenon is 
included in our model, so one can observe the alpha 
structure in the resulting nuclear matter distribution, 
particulary for even-enen nuclei with equal proton 
and neutron numbers. To show this we present in 
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Figure 3: The position and the variance of the wave 
packet for the nucleons in the nucleus ^^O. In the 
first column of the table contained in the figure there 
are the numbers of groups of nucleons, which are cor- 
related in the alpha particle. The second column in- 
dicates the spins and type of nucleons creating given 
pairs of nucleons. The common position and variance 
of the two wave packets is presented. 

The alpha structures obtained in our calculations 
are identical with the structures presented in pa- 
pers [8], [12] up to the silicon and the number of 
alpha-alpha bonds are the same as that proposed by 
[H] (after 160 there is an increment of three bonds 
for each additional alpha particle). Thus we have 15 
bonds for '^^Si instead of 16 8 . The justification 
(but that need not be the final) for this selection can 
be understood from Fig. 7 for alpha-alpha distance 
distribution. 

This distribution shows that the countable alpha- 
alpha bond distances are centered around 3.23 fm up 
to 3.5 fm. In the case of '^^Si, the distances associated 
with 15 bonds lie below this value and the rest has 
distances greater than 4.0 fm. In our model calcula- 
tion, a systematic increase in the binding energy can 
be observed up to ^'^Ni. In this range (2 < Z < 28), 
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Figure 4: The position and the variance of the wave 
packet for the nucleons in the nucleus ^^Mg. Nota- 
tion is the same as in Figure 3. 



Figure 5: The position and the variance of the wave 
packet for the nucleons in the nucleus ^^Ar. Notation 
is the same as in Figure 3. 



systematics follows Bethe's prediction [13], and the 
resulting geometrical arrangement of alpha particle 
has a " balloon structure" (no alpha particles present 
inside the structure) . The model results also indicate 
that the first (lightest) nucleus for which the alpha 
particle appears inside the structure is ^^Zn (see Fig 
8). This results in a departure from the systematic 
increase in energy presented in the Fig 9. 



6 Conclusions 

A new form of the EOS suitable for the MLDM model 
calculations is presented. Preliminary results show 
that the MLDM model is able to reproduce the basic 
properties of atomic nuclei. We see also that MLDM 
fairly well describes the properties of very light nu- 
clei. Because these particles are usually manufac- 
tured with the highest probability this is important 
for a correct description of the reaction dynamics. 

Another important property of the model, result- 
ing from taking into account the mutual interaction 
related to the spin and isospin of nucleons, is the 
emergence of the alpha structures which is particu- 
larly noticeable for even-even nuclei with equal num- 
bers of protons and neutrons. This formation of alpha 



clusters is in line with the considerations on the bind- 
ing energy gain associated with the resulting increase 
in the number of bonds between alpha particles. It 
seems that the model can be helpful in determination 
of the number of bonds. An important prediction of 
the model and strong support for the hypothesis of 
alpha clusters is the appearance of alpha particles in- 
side surface structure for the Z > 28th. This creates 
an additional branch in the scheme of the binding 
energy growth as can be seen in the experimental 
systematics. A detailed analysis of the binding en- 
ergy (as energy per nucleon) for the cases discussed 
above show that there are some additional effects. 
This shows that the number of alpha bonds may dif- 
fer slightly in each individual case. The analysis of 
this behavior as well as restrictions on the applicabil- 
ity of the present simple model will be presented in 
forthcoming paper. 

We need to emphasize the very preliminary char- 
acter of these results. In order to better predict the 
values of the EOS parameters further work is needed 
and we have to perform a global search in which we 
will take into account all available experimental data. 
It would also help to compare the model results, de- 
scribing the basic observables for heavy ion reactions 
(particulary for the compressibility factor Kg) in the 
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Figure 6: The position and the variance of the wave 
packet for the nucleons in the nucleus '^'^Ca. Notation 
is the same as in Figure 3. 



low energy region where the impact of the nucleon- 
nucleon collisions can be neglected. 




bond distance (fm) 



Figure 7: The distribution alpha-alpha intcrparti- 
cle distance calculated with MLDM is plotted (up 
to ^^Ni). Most of the countable distance lies below 
3.5 fm (not a final conclusion). For "^^Si, 16*'* bond 
length is shown. 



0, 77p == only the following terms can appear: 



' = "0 + -■=!;{« + 



5S + ■ 



SSnSn^n + -" 5S Sn Sp'^"^P + -"SSSpSp'^p + ''"iSSi + -'((38^ 
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Appendix A 

In this appendix, it is derived the form of EOS which 
is used in the paper. We start from the expansion of 
the energy density function and we use the symmetry 
which exists in nuclear interactions. According to 
this symmetry, in the expansion (up to 4 order) of the 
e{^, 5,1^,1, rjp) around the point: ^ = 0, (5 = 0, r/„ = 



2/1 1 2 1 2 1 1 

2/1 1 2 1 1 2~\ 

> 1^2''''P'7P + ^='7p1p'7p1p''p + J^5>7p'7p5 + 1''liliVpVp^ ) + 

>7r.^p ( = T7n'7p + =«rj„.7p« + -^iilnVp^^) '^l' 

The symbol e with the index constituting one of 
the variables of the equation denotes a derivative 
of the energy density with respect to this variable. 
We neglect the dependence on (J** and for small spin 
polarization the terms in the brackets containing 
77^, ?7p, Tjnrjp can be neglected in relation to ess, eri„ri„ , 
^VpVp ^ ^^"^ finally one can write: 



1 



eo 



:e«e' 



:ess 



1 ^ 1 X 
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Figure 9: Total binding energy of the nuclei versus 
number of bonds. From ^^C to ^^Ni , it follows 
the same linear relationship, according to surface- 
like structures. From ^'^Zn , slope has a new linear 
relationship reflecting change occuring in the nuclei 
structures, (see text for detail). 



Figure 8: The alpha particle structure of ^"Zn ob- 
tained via MLDM calculation. Noticeable thing is 
that there is an alpha particle inside the balloon like 
structure (the structural change starts from here, see 
the text for detail). 
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2^ivpVp^ + A^iiVpVp^ 



^VnVp \ n^rinVp 



1 



+ T^Hv^np^ + ^Hiv^np^"^ ) (32) 



Replacing symbols of derivates CA.by variables de- 
fined in the following relations: 



i^n 



1 



T8 " 2'«« 



L, 



K, 



CiiO — 2^rjniu ' o ~ 2^^''"''" ' 18 ~ 4^??»M')n 
^yO = 2^'7n')p ' ~o~ ~ 2^iV7zVp ' '7W' ~ '^^iiVnVp 

(33) 
and taking into account the mentioned symmetries of 
nuclear interactions (the derivatives {ejj^rip , ^nnv 



J, 



(e«r,pr)p, e^^„^„) and (e^^^^^^^, e^^^„^„) are pairwise 
equal, this means that the respective symmetry en- 
ergy associated with the spin is identical for protons 
and neutrons), we obtain the equation ([7]). 



Appendix B 

In this appendix we present the nietod of calcula- 
tion of the variance of energy given by the EOS in 
the neighborhood of the center of the packet k. In 
the first step we calculate the average energy density 
of matter at random points with the distribution of 
V(r). 



e(fc) = / e(p,<5,r/„,?7p) •'"/9(r)d^r 



(34) 



If we use the equation of state in the form of an ex- 
pansion ([7]) and (jlip for the isospin part the derived 
term cannot be solved analytically: 



d'^aip ■^p{v)(fr = ai 



Pp 



p-^p{v)d^r 



P 

(35) 
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In order to avoid this problem we can use the fohow- 
ing approximations: 



aj 



f PnPn fe / N >3 ^ ^ f 

J V(r)d-^r ^ai— pn 



V(r)d^r (36) 



where N and A are the total number of neutrons and 
nuclei in the system, respectively. As we can see, 
after removing p from the denominator the integral 
can be solved analytically. Similarly, we can find the 
remaining ingredients (|55)) (with pnPp and PpPp )■ 

In the next step we assume that the 6, r]n, ijp 
are established in the vicinity of the nucleon k as 
S{k), fjn{k), fjp{k) (we replace them by the corre- 
sponding average values) and we determine the ef- 
fective density of Pef{k) based on the equation: 

e{k) ^ e {pef{k),6{k),fin{k),fip{k)) (37) 

In further considerations we replace the set 
Pef{k), 6{k), fjnik), f]p{k) by 

Pptik), Ppi, Pntik), Pniik). To calculate the 
variance <j1{k) we develop the integral expression: 



al{k) = / (g 



e{Pp\,Ppi,p7ii:,Pni)f Mr)d^^ (38) 



around Pp^{k), ppi, Pn\{k), Pnx{k). This allows 
the analytical calculation of the the considered in- 
tegral. In our calculations we limited ourselves 
to the first of order expansion. So, because 
e {Pptik), Ppiik), p„f{k), Pniik)) = e{k) this choice of 
expantion gives a very simple form of the integrand. 
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